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ABSTRACT 

Gravitational microlensing events provide a potentially powerful tool for the study of stellar pop- 
ulations. In particular, they can be used to discover and study a variety of binary systems. A large 
number of binary lenses have already been found by the microlensing surveys and a few of these 
systems show strong evidence of orbital motion on the timescale of the lensing event. We expect that 
more of such binary lenses will be detected in the future. For binaries whose orbital period is com- 
parable to the event duration, the orbital motion can cause the lensing signal to deviate drastically 
from that of a static binary lens. The most striking property of such lightcurves is the presence of 
quasi-periodic features, produced as the source traverses the same regions in the rotating lens plane. 
These repeated features contain information about the orbital period of the lens. If this period can be 
extracted, a lot can be learned about the lensing system even without performing a time-consuming 
detailed lightcurve modeling. However, the relative transverse motion between the source and the 
lens significantly complicates the problem of period extraction. To resolve this difficulty, we present 
a modification to the standard Lomb-Scargle periodogram analysis. We test our method for four rep- 
resentative binary lens systems and demonstrate its efficiency in correctly extracting binary orbital 
periods. 



1. INTRODUCTION 

During the past two decades gravitational microlens- 
ing surveys have transformed from a bold theoretical idea 
( Paczyhski|1986 ) into an important tool for studying the 
stellar population of our Galaxy. Following Paczyhski's 
proposal, four diffe rent mon itoring programs, M ACHO 
( lAlcock et al][l993|), O GLE (lUdalski et al.|1993D , EROS 
flAubourg et al.||l993[ ) and MOA ( |Yock||l998| ), have for 
years collected data on the Galactic Center as well as 
LMC and SMC, detecting many thousands of microlens- 
ing events. Two microlensing surveys that are still cur- 
rently in operation, OGLE-IV and MOA, continue to add 
to this number at a rate of roughly 1000 events per year. 
While the majority of the detected events are consistent 
with being produced by a single-mass lens, a sizable frac- 
tion shows clear evidence of lens binarity ( Alcock et al 



2000| |Jaroszynski et al.||2004} IJaroszynsky et al.||2010l) . 



I'hese are particularly interesting, because m ore co mpli 



cated lightcurves (see e.g. Mao fc Paczyhski||1991[ ) pro 
duced by binary lenses allow us to learn considerably 
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more about these systems (see e.g 
T995||DiStefano &; Pernafl996t . 

The signatures ot lens binarity are most prominent 
when the projected orbital separation between the two 
binary components, a, is of the order of the Einstein ra- 
dius: 



Re 



UGM (D s - Di)Di 
c 2 D„ 



(1) 



where D s and Di are the distances to the source and the 
lens, respectively and M is the total mass of the lens. 
The timescale for the duration of the lensing event is 
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determined by the Einstein radius crossing time 
< E = R E /v, 



(2) 



where v is the relative velocity between the lens and the 
source as viewed by the observer. For solar-mass lenses in 
the Galactic disk is typically on the order of months, 
while the orbital period of the lensing binary, T or b, is 
more likely to be measured in years. As a result, most 
binary lensing events can be modeled while completely 
neglecting orbital motion of the two m asses (see e.g. re- 



cent estimates by Penny et al. 2011a). Nevertheless, a 



handful of microlensing lightcurves do show convin cing 



evidence of binary phase change during the event ( Do- 
| minik|19981|Albrow et al.|2000HAn et al.|2002l|Jaroszyn- 



|ski et al.|2005[|Hwang et al.||2010[|Ryu et al.|2010| . And 

recently there has been an increase of interest in the pos- 
sibility of finding lensing events for which the timescale 
ratio, defined as 

n = t E /T mb , (3) 



is of order of or greater than unity ( Penny et al. 2011b| 
DiStefano fc Esin|[2^TT| ) . 

It the orbital period can be determined for such sys- 
tems, it will provide a new avenue toward extracting the 
lens parameters. Specifically, the period relates the total 
mass to the orbital separation. If the value of a and q 
can also be determined, we have some of the key elements 
needed to specify the binary. In cases when some addi- 
tional information is available, e.g., because the lens is 
located nearby and/or because parallax or finite-source- 
size effects are detected, a full binary solution may be 
obtained. Unlike singl e-lens lightcurves which have a 
simple analytical form ( |Paczyriski| [1986 ) , binary lenses 
can produce very complex magnification patterns which 
cannot be expressed in a closed form. Even for a sta- 
tionary binary lens, the process of fitting an obs erved 
lightcurve can be complex (DiStefano & Mao 1996). The 
addition of orbital motion complicates the fitting process 
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even further. Fortunately, in the limit when 1Z > 1 the 
binarity features observed in the microlensing lightcurves 
show very strong periodicity which allow T or b to be deter- 
mi ned via simple Lomb-Scar gle (LS) periodogram analy- 
sis ( DiStefano fc Esm||2011 ). In this paper, however, we 
are primarily interested in the regime when TZ < 1. For 
these longer-period binaries the interplay between the 
orbital motion and the projected motion of the source 
causes the period det ermined via LS analysis to differ 
significantly from T or b (Di Stefano & Esin 2011 ). Here we 
describe a modified timing analysis method which allows 
us to compensate for this relative motion and extract the 
correct binary period of the lens. 

2. MICROLENSING BY A ROTATING BINARY 
LENS 

A microlensing event produced by a single lens can be 
completely characterized by three parameters: -Re, the 
Einstein radius of the system; b, the distance of clos- 
est approach between the lens and the projection of the 
source onto the lens plane (measured in units of -Re); 
and v, the transverse velocity of the source with respect 
to the lens. In general, we can express the separation 
between the lens and the source projection (measured 
again in units of -Re) as 



u(t) = Vb 2 + l(t-t )/t E }< 



(4) 



where to is the time of closest approach. The result- 
ing lightcurve is described in terms of the amplification 
function, defined as the ratio of the observed lensed and 
unlensed fluxes, which takes the form 



u\J u 2 + 4 



(5) 



When the lens is a binary, in addition to the total 
mass of the system, M = Mi + M2, we need to specify 
the mass ratio q = M\jMi and the dimensionless semi- 
major axis of the binary, a — a/ Re, to fully describe 
the microlensing event. For simplicity, in this paper we 
consider only face-on, circular binary orbits, so a is equal 
to the binary separation at all times. Since the lens is 
rotating, two more parameters are necessary: the value 
of the binary phase angle (shown as angle <fi in Fig. [I]) 
at the moment of closest approach, </>o, and the direction 
of binary rotation with respect to the relative motion 
vector. For the source velocity shown in Fig. [II we define 
a binary to be prograde if it is rotating clockwise and 
retrograde otherwise. 

Since general relativity is a non-linear theory, the light 
curve produced by a binary lens is more complicated than 
a simple superposition of two amplification functions and 
cannot be described by an analytical formula. However, 
we can calculate the amplification value nume rically as a 



function of so urce p osition in the lens plane ( [Schneider 
& Weiss||1986||Witt||1990|). The result depends only on 



a and q. Given b and v, we can calculate the trajectory 
of the source projection in the plane of binary lens and 
construct the resulting lightcurve. 

Fig. pla) shows a sample lightcurve produced by our 
standard binary lens system (described in detail in <Q. 
Over-plotted in dashed red line is the best-fit single-lens 
light curve (Eq. [5]) . Overall shapes of the two lightcurves 
are very similar, which is not surprising for a binary with 




Figure 1. Geometry of the lensing system. The black dot shows 
the projection of the source location onto the lens plane. The 
trajectory of the source in the lens plane is given by the dashed 
line. Note that we refer to the vertical line connecting the center 
of mass of the binary lens and the point of closest approach of the 
source as a CoM-b-axis. 

a fairly low value of a = 0.25. However, subtracting the 
single-lens fit from the binary lightcurve reveals a clear 
oscillatory pattern in the residuals, shown in Fig. [2jb) . 




400 600 800 
t (days) 



1000 1200 



Figure 2. (a) The black solid line shows a sample theoretical 
light curve with b = 1.0 produced by our standard binary lens- 
ing system (see ijljl rotating in retrograde direction. The dashed 
red line shows the least-x 2 fit single-lens lightcurve Agi^(t). (b) 
The residual lightcurve, defined as <5(t) = A(t) — AslMi clearly 
shows the effects of lens binarity. The oscillatory pattern contains 
information about the orbital motion of the lens. 

To understand the relationship between this periodic 
pattern in the residuals and the orbital motion of the 
lens, it is instructive to examine the trajectory of the 
source in the lensing plane. In Fig. [3] we show this 
trajectory superimposed onto the contour plot of the 
absolute value of the residual magnification \S(x, y)\ — 
\A(x, y) — ^4sl(^, y)\- Here x and y are coordinates in the 
lens frame, scaled to -Re- The residual contours exhibit a 
nearly-symmetrical four-lobe shape. Note that the longer 
lobes (lined up with the y-axis) correspond to positive 
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5(x,y) values , while the shorter ones have 5{x,y) < 0. 
For face-on circular orbits a is time-independent and the 
amplification map is static. In the co-moving and co- 
rotating frame of the lens, the source appears to spiral 
in as it moves towards the point of closest approach and 
spiral out as it moves away from the lens. If the source 
were stationary, one full rotation of the binary lens would 
result in two full oscillations of the residual lightcurve. 
Thus, the periodicity of the residual signal is related to 
the orbital period of the binary lens by T ~ T or b/2. Note 
that the relation is not exact due to the modification 
caused by relative transverse motion between the source 
and the lens. 



where ui n 
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Figure 3. Both panels show the contour plot of the absolute 
residual magnification function, = |A(:e,jj) — A^{x, y)\, 

produced by the standard system (described in 941. The contour 
map is centered on the center of mass of the binary lens and the 
axes are in units of ijg . The two binary components are located on 
the ce-axis at x = ±a/2 = ±0.125. The horizontal elongated petals 
correspond to negative residuals and vertical petals correspond to 
positive residuals. The legend in panel (a) specifies the minimum 
values of \&(x, y)\ that correspond to different contour colors. As 
the source moves past the lens, its trajectory in the lens plane 
spirals in and out, producing the oscillatory residual pattern seen 
in the lower panels. Note that the colors of the trajectories and 
the residual curves are matched. Panel (a) shows the trajectory of 
the source with b = 1.0 for the prograde binary lens, while panel 
(b) corresponds to the retrograde case (the lightcurve is identical 
to the one shown in Fig. [2j. 



To get a quantitative understanding of the interplay 
between source motion and binary rotation, we need to 
derive an equation for the time evolution of angle 9, as 
shown in Fig. [T] Assuming that the source moves with 
velocity v in trie lens plane, its distance to the CoM-b- 
axis, scaled to i?E, takes the form 



t-tg 



(6) 



The angle from the CoM-b-axis to the source is then 
given by 

, . w t — to , . 

7m = — arctan — = — arctan . (7) 

o toe 

Notice that 7(f) > when t < to and j(t) < when 
t > to- The phase angle from b-CoM-axis to the axis of 
the binary system is 



^orb = i s defined as the orbital frequency 
and "— " sign is adopted for prograde motion while "+" 
sign corresponds to retrograde motion. Finally, the angle 
between the source position and the axis of the binary is 
then 



0(t) = tj>(t) - 7 (i) 



6q =F ^o-cht + arctan 



t-tp 
bt v . 



(9) 



<t>{t) = </>o T u orb t, 



(8) 



The distance from the source to the center of mass in 
the lens plane, u(t) is given by Eq. |4j). Note that the 
position of the source in the co-rotating lens plane is 
simply 

x(t) = -u(t)cos9(t), y(t) = u(t)sw.9(t). (10) 

A natural consequence of Eq. ([9| is that the angu- 
lar velocity of the source in the lens plane, 9, is smaller 
for a prograde system than that for a retrograde sys- 
tem. This effect manifests itself clearly in the apparent 
period of residual signals shown in Fig. |3j which com- 
pares the trajectories of the source for a prograde and 
retrograde system with otherwise identical physical pa- 
rameters. The retrograde residual lightcurve shows much 
faster oscillations due to having a larger value of 9. 

3. TIMING ANALYSIS 

We now address the question of how to extract the in- 
formation about the orbital periodicity of the binary lens 
from the observed microlensing lightcurves. We start by 
generating a theoretical binary lightcurve, add noise to 
simulate real survey data, subtract the dominant single- 
lens signal, perform spectral analysis on the the residuals 
and examine the accuracy of the inferred orbital period. 

3.1. Modeling the Observations 

We can generate a theoretical light curve of arbi- 
trary duration. To approximate realistic observing con- 
ditions, we take the photometric uncertainty to be 1% 
for most of our calculations. (We discuss the results for 
high-precision photometry in <^6j) At u = 3 Eq. ([5| 
gives Asl ~ 1.01, so the overall amplification level is 
marginally distinguishable from random noise fluctua- 
tion. Based on this, we define the span of each lensing 
event to be from 3£e before the time of closest approach 
to 3^e after the closest approach. For all our lightcurves 
we assume a fairly conservative average sampling interval 
of At = 2 days. With this average rate, we sample the 
theoretical light curve at random times {tj} and add ran- 
dom Gaussian noise with mean and standard deviation 
0.01 to each data point. 

Finally, we fit Eq. (15) to our simulated lightcurve, both 
to obtain the residuaTlightcurve S(t) as well as to deter- 
mine our best guess values for b,to,tE, which are neces- 
sary for the subsequent timing analysis. We found that 
the fitting procedure works much better after smooth- 
ing out the sharp features, if any, produced by caustic 
crossings. 

3.2. Modified Lomb-Scargle Analysis 

Once we have the residual lightcurve we are ready to 
perform spectral analysis to extract the underlying pe- 
riodicity of the binary. We base our period extraction 
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method o n the classical LS periodogram analysis (Scar- 
gle 1982), a modification of the Fourier transform de- 
signed tor unevenly sampled data that has the advantage 
of time-translation invariance. 

For time series {S(tj)} with N data points, the LS 
periodogram as a function of frequency uj is defined as 



where Cj and Sj are given by 
Cj = cos(cu (tj 
Sj = shi(ui(tj - 



■r)), 
r)), 



(11) 

(12) 
(13) 



and r is defined as 



r = — tan 
2u 



N 

sin 2uit j 

3=1 

N 

cos 2uitj 

3=1 



(14) 



Unsurprisingly, the unmodified LS periodogram anal- 
ysis does not in general yield the correct binary period, 
since the frequency of the observed oscillatory signal, oj, 
is affected by the relative source- lens motion (see Fig. |3| 
and is therefore not equal to the the binary orbital fre- 
quency, Worb- Figure EJjb) shows the LS power-spectrum 
for a microlensing ligntcurve with b = 0.4 generated by 
our standard lensing system. The peak power corre- 
sponds to a period of 266 days, very different from both 
the true orbital period of the system, T or b = 372.5 days, 
as well as T or b/2 (which we expect to come closer to 
reflecting the true variability timescale of the lensing sig- 
nal due to the bi-fold symmetry of the magnification map 
produced by this binary, as illustrated by Fig. [3]). 

In <|2] we argued that the progress of the lensing event 
is governed by the time evolution of angle 8, defined by 
Eq. ([9]), which takes into account both binary rotation 
and source motion. Thus, we propose to modify the os- 
cillating coefficients Cj and Sj, essentially replacing uitj 
with 0:1, i: 

Cj = cos[2{w orb (i, - r) - j(tj)}}, (15) 

Sj = sin[2{o; orb (t j - r) - y(tj)}]. (16) 

We dropped the <po term in Eq. (J9| since it simply pro- 
duces a time translation of the signal without affecting 
the results of the LS analysis. The extra factor of 2 origi- 
nates from the fact that the binary magnification pattern 
repeats itself twice over one binary revolution, so the pe- 
riod directly observed in the data would be roughly half 
of the actual binary period, as we already pointed out 
above. Note that the definition of t remains the same 
apart from replacing ui with u; or b in Eq. ( |14[ ) . Finally, the 
approximate values of b, to and £e necessary for evaluat- 
ing 7 were obtained when fitting Eq. ([5]) to the original 
lightcurve. 

Unlike the standard LS periodogram analysis for which 
only positive ui values have meaning, u; orb can have ei- 
ther sign. Following the sign convention in Eq. ([9]), when 




Figure 4. Periodogram of a lightcurve with b = 0.4 produced by 
the standard system in retrograde rotation. Panel (a) shows the 
results produced by the modified LS analysis. The highest power 
peak corresponds to the period of 362 days, very close to the actual 
orbital period of the binary shown with a dashed line. Panel (b) 
displays the results of the standard LS periodogram analysis (here 
the power spectrum for negative period values is simply a mirror 
image of the positive part) . The best fit period of 266 days is clearly 
very different from either T OT \^ or T or t,/2 shown by dashed lines. 



interpreting the results of the modified periodogram 
w or b < corresponds to prograde systems while w rb > 
corresponds to retrograde systems. 

Fig. [3] illustrates the difference between the results of 
the standard (panel [b]) and modified (panel [a]) peri- 
odogram analyses. Our new method produces a sharp 
power peak corresponding to the period T = 27r/w or b = 
362 days, within 3% of the actual binary period. In addi- 
tion, a positive period value correctly identifies the lens 
as a retrograde binary. 

3.3. False Alarm Probability 

The final issue we need to address is the how to de- 
termine the minimum power level, P C ut, such that any 
peak with a power above P cu t is unlikely to be produced 
by chance. Scarglc (1982 ) derived a formula relating the 
false alarm probability p{ and P cu t: 



Pcut = -In[l-(1-P/) 1/J H (17) 

where N u is the number of independent frequen cies 
searched. For pf <C 1, we can approximate Eq. (JT7J) 
as 

P cut » HN^/Pf) = ln(i\y + ln(l/ P/ ). (18) 

For a given value of pf, P cut depends only on N u . There 
is no clear consensus in the literature on how to calculate 
Nu> for a given set of data, sinc e the precise answer de- 
pends on the sampling method ( Gilliland fc Fisher|1985 



Baliunas et. al.|1985 Home & Baliunas 1986). However 
it is clear that iV^ must be on the order ot the number 
of data points, N. We can thus rewrite Eq. (18) as 



P cut - ln(N) = ]n(N u /N) + m(l/ P/ ), (19) 

and determine the term hiN^/N empirically, using the 
same sampling scheme as we do for our light curves. 
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In order to do this, we simulated 5000 time series, ran- 
domly sampling pure Gaussian noise with mean and 
a = 1. Each time series consisted of N data points, with 
N randomly chosen between 10 and 1000 using a uniform 
distribution in log space. We then applied the standard 
LS analysis to each time series by searching through fre- 
quencies from 2tt to Nir, and examined the distribution 
of Pmax — ln(7V) values, where P max is the power corre- 
sponding to the highest peak in each time series. We 
found that P max — ln(iV) tended to increase as we in- 
creased N u but above N u = 16N the distributions re- 
mained virtually unchanged. This worst case scenario is 
shown in Fig. [5] Based on this simulation, if we choose 
-Pcut = 6 + ln(7v) as our significance criterion, only 18 
out of 5000 time series register P max above this thresh- 
old. This corresponds to the false alarm probability of 
18/5000 = 0.0036, which compares well to p f = 0.0025 
we get from Eq. (181 if we set N u = N. We adopted 
F cut = In (AT) + 6 as our detection threshold, i.e. only 
the peaks with power exceeding P cut were treated as real 
detections. 




Figure 5. Distribution of P max — In(jV) values for 5000 Gaussian 
random noise time series. The probability of detecting a spurious 
peak with a power P > 6 + ln(7V) is less than 0.36%. 



4. THE STANDARD SYSTEM 

Having established our method for timing analysis, we 
first test it for what we call our standard binary lensing 
system. It consists of two 0.5M Q stars separated by a 
distance a — 1 AU. We would expect this type of binary 
system to be fairly common since O.5M stars lie near 
the peak of the IMF. We take the source to be near the 
Galactic Center at a distance of 8kpc and set the lens 
at 4kpc. For the relative transverse velocity we adopt 
the value v = 60km/s. For these parameters, the ratio 
of the binary separation to the Einstein radius is a = 
0.25 and the timescale ratio is 1Z = 0.31. These are 
fairly conservative choices. On the one hand, a is low 
enough that the signatures of binarity will not be very 
prominent, as d emonstrated in Figs. [2] and [3] (see also 

i|20fT|. In addition, a 



the discussion in Di Stefano fc Esin|: 
relatively low value of the timescale ratio puts this binary 
firmly in the regime where both the orbital motion and 
the transverse motion of the source will play important 
roles in the formation of the lightcurve, making it an 
ideal testing system for our timing analysis. 

Our goal is to investigate the success of our proposed 
period extraction method for this system as a function 



of the impact parameter b. We want to determine how 
weak must be the overall amplification of a microlens- 
ing event when the periodicity in the lightcurve can no 
longer be detected. To this effect, we vary the value of 
b from 0.05 to 4.0. Since low b means stronger events 
and potentially more identified lightcurves in the exist- 
ing lensing archives, we adopt a finer mesh for b < 0.8 
and a coarser mesh for larger values of b. 

For each value of b we generated 2000 lightcurves 
equally split between prograde and retrograde rotation 
directions. For each lightcurve, the initial phase angle 
0o was randomly chosen from the interval [0, 2tt]. Each 
lightcurve was then analyzed using our method to deter- 
mine (1) whether the power spectrum contains at least 
one peak with P > P cu t, as defined in the previous sec- 
tion; and (2) whether the period corresponding to the 
peak with the highest power falls within 10% of the true 
orbital period of the lensing binary. 

4.1. Period Detection Rates 

The results are summarized in Fig. [6] As expected, 
the overall period detection rate (shown as a dotted black 
line) falls off with increasing b. This makes perfect sense 
because as the lensing event becomes weaker, the binarity 
signal becomes buried in the noise. Nevertheless, it is 
significant that most lightcurves are identified as periodic 
past 6 = 3. 




Figure 6. The efficiency of our timing analysis for the standard 
system. The dashed black line shows the period detection rate, i.e. 
the fraction of the simulated lightcurves for which of a significant 
periodic signal is detected. The solid red line shows the correct 
period detection rate, i.e. the probability of extracting a significant 
orbital period that is within f0% of the true period of 372.5 days. 



Of more interest to us is the correct period detection 
rate (shown as a solid red line), i.e. the probability of 
extracting a significant orbital period that is reasonably 
close to the true period of the binary. For the purposes of 
this paper we defined a detected period to be "reasonably 
close" if it lies within 10% of the true binary period of 
the lensing system (this also includes having the correct 
sign, i.e. distinguishing between prograde and retrograde 
rotation). This rate is falling faster with increasing b 
than the overall detection rate, but even at b = 3, we 
are still correctly identifying periods for 30% of all the 
lightcurves. Figure [7] shows the distribution of detected 
periods for lightcurves with b — 1.0. It is clear that for 
most of these events we can use our method to determine 
the binary period to an accuracy of a few percent. 

While the behavior of the detection rates for b > 1.0 is 
easily understood, there are two peculiar features evident 
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Figure 7. Distribution of periods detected for b = 1.0. 



for the lightcurves with b < 1.0. Firstly, for b < 0.2 there 
is a steep drop in both detection rates. In addition, in the 
region 0, 2 < b < 1.0 there is a ~ 50% drop in the correct 
detection rate, while the overall detection rate remains 
near 100%. We investigate these in the remainder of this 
section. 



analysis and since it is never repeated no significant pe- 
riod is detected in the data. 

A simple way of dealing with this issue is to simply 
remove the high-frequency part of the lightcurve so that 
the global low-frequency signal in given more weight in 
the timing analysis. To determine which part of the 
lightcurve to remove, we estimated the average radius 
of the red region to be ~ 0.2i?E, and removed the cor- 
responding portion of the lightcurve, i.e. the region 
with the overall magnification higher than Asl(0.2) sw 5, 
where is given by Eq. ^ . The detection rates with 
and without this modification to the timing analysis are 
shown in Fig. [9] Cutting out high magnification points 
appears to be highly effective; it restores the period de- 
tection rate to near 100% and raises the correct detection 
rate to at worst 50% (Fig. [9la)). Our investigation of 
other systems (£[5]) showed that his high-frequency sig- 
nal contamination during close-approach events appears 
to be an universal problem. Thus, any results we show 
from this point on include this step in the timing analy- 
sis. 



4.2. High Frequency Signal in Low-b Events 

To diagnose the precipitous drop in the detection rates 
for b < 0.2, we examine the trajectory and the cor- 
responding residual lightcurve curves of an event with 
b = 0.1 (see Fig. M. 
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Figure 8. The source trajectory and the corresponding residual 
lightcurve for a prograde standard system with b = 0.1. 
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Figure 9. Comparison of the detection rates with (panel [a]) and 
without (panel [b]) removing the high magnification (A > 5) por- 
tion of the lightcurves. The dotted black and solid red curves are 
defined as in Fig. [6] 



4.3. Prograde-Retrograde Period Confusion 

To investigate the drop in the correct detection rate 
in the region 0.2 < b < 1.0, we examined the statistics 
for prograde and retrograde systems separately. Fig. [10] 
demonstrates that this effect is mainly due to the pro- 
grade systems. 




Figure 10. The detection efficiency for the standard system for 
lightcurves with prograde (panel [a]) and retrograde (panel [b]) 
rotation. 



The problem arises because the overall shape of the 
residual curve shows a single strong high-frequency oscil- 
lation near the peak of the microlensing event. Relating 
the residual curve to the trajectory, we can see that it 
is due to the crossing of the red region between the two 
binary companions. Because of its relatively large ampli- 
tude, this high-frequency feature dominates our timing 



Next we examined the statistically significant though 
incorrect periods detected for the prograde lightcurves. 
It turns out that for many of these events the highest 
peak in the periodograms corresponds to a longer retro- 
grade period ~ +500 days instead of the correct prograde 
period of —372.5 days, The correct period is also de- 
tected, but at a slightly lower power than the retrograde 
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Figure 11. Periodogram of a prograde standard system with b = 
0.4. The highest peak corresponds to a longer retrograde period 
+491 days while the true prograde period —372.5 days has a lower 
power. The dashed line shows the location of the true system 
period. 



period. A per iod ogram of such a prograde lightcurve is 
shown in Fig. [TT] 

This confusion originates when subtracting the fitted 
single-lens lightcurve. Even a slight discrepancy between 
the actual and fitted values for 6, £e and to introduces 
an extra bump near the peak of the lightcurve. This ex- 
tra feature causes the prograde residual lightcurves to be 
mistaken for a retrograde one, since the introduction of 
the extra peak in the prograde case will mimic the effect 
of higher frequency oscillat ions characteristic of retro- 
grade lightcurves. In Fig. [12] we plot the probability 
that one of the two highest-power peaks lies within 10% 
of the correct period (dashed blue line), in addition to 
the two detection rates we have been discussing so far. 
This new detection rate shows that indeed most of the 
missing periods for low-6 lightcurves appear in the tim- 
ing analysis as second highest power peaks. Interestingly, 
the detection rate of periods for higher b lightcurves also 
significantly improves is we are willing to settle for two 
possible answers for the binary period. 



5. OTHER EXAMPLE SYSTEMS 

We now examine the efficiency of our timing analysis 
technique for four other example lensing systems. For 
all of them we take Di and D s to be the same as for 
our standard binary. The rest of their physical parame- 
ters are summarized in Table [l] an d th e resulting period 
detection rates are plotted in Fig. [13] Note that panel 
(a) shows the results for our standard system for ease of 
comparison. 



system 


M/Mq 


g 


T orb (d) 


v (km/s) 


a 


n 


Standard 


1 


l 


372.54 


60 


0.25 


0.3135 


Unequal 


1 


0.5 


372.54 


60 


0.25 


0.3135 


Tight 


1 


1 


172.76 


60 


0.15 


0.676 


Wide 


1 


1 


752.31 


60 


0.4 


0.155 



Table 1 

Summary of relevant parameters for the four example systems. 



1 2 3 1 2 3 4 




1 2 3 1 2 3 4 

b b 




Figure 12. Detection efficiency in standard system. The blue 
line represents the probability of either one of the highest peak 
and second highest peak lie in within 10% of the correct period of 
the system. 



It turns out that this confusion between prograde and 
retrograde rotation is specific to our standard system 
rather than universal to binary lensing systems. In most 
of the other example systems that we discuss in ^\ this 
issue hardly arises except for one that is most similar to 
the standard binary. 



Figure 13. Period detection rates for four systems summarized 
in Table 1. 



5.1. Unequal Mass Binary 

This binary system is essentially equivalent to the stan- 
dard binary lens system we discussed in the previous sec- 
tion except the masses of the two companions are |M© 

and §M©, respectively. The effect of having q ^ 1 is the 
change in the shape of the residual amplification map; 
as illustrated in Fig. [14J the spatial asymmetry of mass 
distribution about the y-axis leads to the tilting of the 
"petals" in the residual pattern. As the source follows 
its spiraling trajectory in the lens plane, this asymmetry 
will certainly affect the spacing of the peaks, as illus- 
trated in the lightcurves shown in the lower panels in 
Fig. [l4j and can complicate the process of extracting 
the orbital period. However, it is encouraging t hat the 
detection efficiency for this system, shown in Fig. 



13 >), 



is only marginally diff erent than for the standard equal- 
mass binary (Fig. [L3[a]). In fact, the overall shapes of 
the detection rates for the two cases are essentially iden- 
tical when b < 1. At larger values of b, for which the 
asymmetry is more pronounced, the unequal-mass sys- 
tem shows slightly lower detection rates, but the drop 
does not exceed 5-10%. 
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Figure 14. The magnification map for the unequal-mass binary 
system. Since the two masses are not the same, the magnification 
map is not symmetrical with respect to the y-axis. Also plotted 
are the source trajectories with b = 1.0 for both prograde (panel 
a) and retrograde (panel b) binary rotation directions, as well as 
the corresponding residual lightcurves. 



5.2. Tight Binary 

This system differs from the standard case only in its 
smaller binary separation. With a — 0.15, the charac- 
teristic binarity features in the residual lightcurves hav e 
very low amplitudes (see e.g. DiStefano & Esin 2011). 
On the other hand, smaller binary separation also means 
smaller orbital period. As a result, the timescale ratio is 
more than twice as large for this system as for our stan- 
dard case, so the periodicity is easier to detect even in 
shorter duration lightcurves 

it 



Looking at the detection rates shown in Fig. 13 ( 
is clear that our results are dominated by the decrease 
in a, the parameter which determines the amplitude of 
the features in the residual lightcurves. At best, we can 
detect periodicity in ~ 60% of all lightcurves and the 
probability of detection drops nearly to zero for b > 1.5. 
However, because of the relatively high value of TZ, when 
a period is detected the chances of it being correct are 
very high. The only exceptions are lightcurves with b < 
0.1, where the correct detection rate drops to about half 
of the total rate, just like it does for the standard system. 
We show in ^6]that the detection rates for this system are 
very sensitive to the photometric precision of the data; 
with sma ller noise the detectability can become very high 
(seeFig.pt). 



correct detection rate drops significantly, reaching a max- 
imum of 60% near b = 0.5 and dropping off to zero for 
6>3. 

If we investigate the prograde and retrograde 
lightcurves separately, we can see that the gap between 
the two detection rates is caused predominantly by the 
prograde systems. In fact, taken by themselves, prograde 
lightcurves have the correct detection rate close to zero. 
The reason for this is simple. At TZ = 0.155 this binary 
straddles the boundary of detectability since its orbital 
period is slightly longer than the length of the microlens- 
ing event. In prograde binaries the effective period of 
the residuals is larger than T or b, so detection becomes 
virtually impossible. For retrograde binaries, the effec- 
tive period is smaller than T or b and can therefore be still 
detected. 

6. LOW NOISE DETECTION RATES 

While monitoring observations from the ground are un- 
likely to have photometric errors smaller than ~ 1%, 
space-based missions, such as Kepler, can do roughly 
1000 times better. It is highly likely that multiple mi- 
crolensing events will occur in the Kepler field over the 
lifetime of the mission, so it is worthwhile to consider 
how our period detection method would fare for data 
with photometric precision on the order of 0.001%. To 
that end, we have computed the period detection rates 
for our four systems in this low noise limit. Since high 
precision photometry would allow us to detect much 
smaller deviations, we expanded the length of our simu- 
lated lightcurves to 14£e, instead of 6^e we adopted for 
higher noise levels. 
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5.3. Wide Binary 

We now consider the effect of having a larger binary 
separation while keeping all the other parameters the 
same. With a = 0.4, this lensing system will produce 
lightcurves characterized by much more pronounced de- 
viations from the single-lens form. However, larger a im- 
plies longer T or b and therefore smaller TZ, making period 
detection more difficult. 

The results for the detection rates are shown in Fig. 
13 'd). In contrast to the small-separation system dis- 
cussed in the section above, the overall detection rate 
for the wide binary lightcurves shows only a moderate 
decrease compared to the standard case. However, the 



Figure 15. Period detection rates for our four systems with 
0.001% photometric uncertainties. 

Fig. [15] shows our new detection rates for the four 
systems. It is clear that the main effect of increasing 
photometric precision is a significant increase in detec- 
tion rates, especially at high b. This result makes perfect 
sense: at large b, the overall lightcurve amplification as 
well as deviations from the point-lens form become small 
enough as to be non-detectable from the ground, but 
perfectly discernible from space. 

Note that the overall rate of period detection now re- 
mains flat at essentially 100% all the way beyond 6 = 4. 
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Examining retrograde and prograde lightcurves sepa- 
rately, we again find that most of the gap between the 
overall and correct detection rates is attributable to pro- 
grade rotation. 

The systems affected the most by the change in noise 
levels are the tight and wide binaries. The results for the 
former are particularly striking; where before we could 
detect at most 60% of the periods, we now can detect 
the correct period for any lightcurve. The main reason 
for this improvement is again entirely due to the fact 
that very small amplitude signal can now be clearly de- 
tected. The value of a for this system places the binary 
features in the microlensing lightcurves at the edge of 
detectability with 1% errorbars, but this is not the case 
with improved precision. 

For the wide binary, the periods for prograde 
lightcurves were largely not detectable which higher noise 
levels since the length of the event was simply too long 
compared with the orbital period of the lens. Improved 
sensitivity allows us to follow the event for a much longer 
time, effectively decreasing the minimum value of TZ for 
which the periods can be reliably detected. 

While the minimum value of TZ for which our timing 
analysis can be applied is set by the photometric preci- 
sion of the data, we can place a firmer limit on the upper 
value of TZ for which our method is important. Clearly, 
when TZ » 1, the lens-source relative motion becomes 
irrelevant and standard LS analysis should yield the cor- 
rect per iod (up to a factor of 2) for such a binary lens. 
Fig. fl6] illustrates when the standard method results be- 
gin to deviate significantly from the actual orbital period. 
The black and red curves show the timing analysis results 
for a series of lightcurves, all with 6 = 1, produced by 
lensing binaries with 1Z spanning the range between 0.1 
and 2. All the other parameters were chosen to be the 
same as for our standard system with the exception of a, 
which w as a djusted to vary T or b and therefore 1Z. 

Fig. [16] shows that our proposed timing analysis 
method is very successful in extracting the correct orbital 
periods for the entire range of timescale ratios above the 
minimum set by the length of the lightcurves^ while the 
standard LS periodogram algorithm clearly begins to fail 
(i.e. the detected period differs from the actual value by 
more than 10%) when TZ < 1.0 for prograde lenses and 
1Z < 0.8 for retrograde lenses. 

7. DISCUSSION AND CONCLUSIONS 

When the Einstein crossing time is comparable to the 
orbital period of a binary lens, microlensing lightcurves 
display quasi-periodic features which can be used to de- 
termine T or b- In the limit when TZ > 1, the standard 
LS periodogram analysis can be used to find the pe- 
riod. However, in this paper we focus on the regime 
with TZ < 1 when the variability timescale in observed 
lightcurves is set by a combination of binary rotation 
and relative source-lens motion. We have shown that 
this regime includes specific systems with orbital periods 
in the range from months to years. For such lightcurves 

5 Note that for high values of TZ, a becomes very small, but the 
periods are still detectable with the low-noise data. This would not 
be the case for 1% photometry; however, other lensing parameters 
can be adjusted to e nsure that the value of a remains within a 
detectable range (see |Di Stefano fc Esin|20iT| ) . So our conclusions 
for TZ are still applicable. 




Figure 16. The ratio of the detected to actual orbital periods 
of the binary lens as a function of TZ for LS (black triangles) and 
modified LS (red squares) periodogram methods. The results for 
the former were divided by a factor of 2 to compensate for the sym- 
metry of the amplification pattern. Top and bottom panels show 
the results for retrograde and prograde lightcurves, respectively. 



the standard periodogram yields incorrect periods (see 
Fig. |4j). Here we proposed a modification to the timing 
analysis method designed to compensate for the source 
motion and extract the correct orbital period. 

We tested our new timing analysis method on simu- 
lated lightcurves for four different binary lens systems 
and calculated period detection rates for a wide range of 
impact parameter values. The results are very encour- 
aging. As long as the orbital period is detectable (i.e. 
not longer than the length of the simulated lightcurve) 
our proposed method finds and identifies it correctly in 
a large fraction of cases. 

We find that the period detection rates are determined 
primarily by the photometric precision and three interde- 
pendent parameters, a, TZ and b. Perhaps unsurprisingly, 
the photometric precision is the most important. It de- 
termines the minimu m value of a for which b inarity of 
the lens is detectable (|Di Stefano k. Esin||2011 ) and sets 
the maximum value of b tor which a microlensing event 
is distinguishable from the noise. It also essentially sets 
the maximum practical length for a lightcurve which in 
turn determines the maximum detectable orbital period 
and therefore places a lower limit on TZ. 

We consider two levels of photometric uncertainty: 1%, 
characteristic of ground-based data, and 0.001%, relevant 
for dedicated space missions such as Kepler. In the first 
case, the reasonable event duration is 6£e (beyond this, 
the wings disappear into the noise) and, correspondingly, 
the detection rate drops off for TZ < 1/6, as illustrated in 
Fi g. |16| and well as by our results for the wide binary lens 
in §5.3| High-precision photometry allows us to signifi- 
cantly sxpand the parameter space for which the orbital 
periods can be reliably detected with our method. 

The results we discuss in this paper are based on the 
analysis of lightcurves produced by binary lenses in face- 
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on circular orbits. We are now in the process of extending 
our calculations to include elliptical and inclined orbits. 
In these more complex situations, the projected binary 
separation, and with it the shape of the lensing magnifi- 
cation pattern, will vary with time during the microlens- 
ing event. This effect will certainly affect the spacing of 
the features in the residual lightcurves, possibly affect- 
ing the extraction of the orbital period. Our preliminary 
results indicate that the orbital period should still be de- 
tectable in a large fraction of lightcurves, but the full 
report will be the subject of a follow-up paper. 

Finally, we want to point out that in its present form, 
our timing analysis method works best for binary lenses 
with q not very different from unity. It is predicated on 
the assumption that the oscillatory features repeat twice 
per orbit. The fact that period detectability does not 

al and [b]) 
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change much from (7 = 1.0toq = 0.5 (Fig 
shows that some asymmetry in the magnification pattern 
can be easily tolerated. We tested lower values of q for 
the unequal-mass system and found that the detectabil- 
ity decreases by a factor ~ 2 for q = 0.25 and drops 
below 20% for q = 0.1. It is clear that in the plane tary 
regime (i.e., when q << 1) a factor of 2 in Eqs. (151 



and ( 16 1 does not apply, since the magnifi cation pat- 
tern is no longer even remotely symmetrical ( Di Stefano 



2011 1. Fortunately, such systems produce visibly differ 
ent light curve morphologies characte rized by very spiky 



rathe r than sinusoidal features (see Di Stefano & Esin 
|2011[ for some examples of lightcurves produced by low- 
q binary lenses), and so can be flagged. We are now 
working on ways to extend our method to such \ow-q 
systems. 
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